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The process e"*"e~ —^ ZH allows to measure the Higgs boson in the recoil mass spectrum 
against the Z boson without any assumptions on the Higgs boson decay. We performed 
a full simulation and reconstruction of e^e~ — > ZH using the Mokka and Marlin 
packages describing the LDC detector. The Z is reconstructed from its decays into 
electrons and muons. The mass of the Higgs boson is set to 120 GeV. Assuming a 
centre-of-mass energy of 250 GeV and an integrated luminosity of 50fb~^ the Higgs 
boson mass and the Higgs-strahlung cross section can be measured with a precision of 
120 MeV and 9%, respectively. 

1 Introduction 

The Standard Model (SM) [1] predicts one Higgs boson as a remnant from the spontaneous 
symmetry breaking mechanism [5] . This mechanism allows fermions and the W and Z bosons 
to acquire their masses by interaction with the Higgs field. To discover the Higgs boson, 
therefore, is of crucial interest to complete the SM. Electroweak precision measurements 
suggest the mass of the Higgs boson to be of the order 0(100 GeV). Direct searches at 
LEP have set a lower mass limit of 114 GeV. A Higgs boson with a mass above 114 GeV 
will be accessible in the experiments at the Large Hadron Collider (LHC). However, to be 
sure about the nature of the particle found, it is necessary to measure its properties such as 
mass, width, charge, spin, parity, couplings to other particles, and self-couplings to test the 
internal consistency of the SM, or to find hints for new physics. 

For the determination of at least some of these quantities, LHC will not be sufficient. 
At the future International Linear Collider (ILC), we will have the chance to investigate 
the properties of new particles with high precision in all details. This e'^e~ collider with a 
center-of-mass energy up to 1 TeV provides a well known initial state, a very clear signature 
for events of e"'"e" ZH, and due to its high luminosity sufficient statistics for precision 
measurements . 

In the Higgs-strahlung process e+e~ ZH, we can investigate the coupling of the Higgs 
boson to the Z boson and determine the Higgs boson mass with high precision in a relatively 
model- independent way using the recoil mass spectrum against the Z, 

"^recoil = S + ^di-lopton ~ 2 • i?di-lopton ' \/s , (1) 

where s is the square of the centre-of-mass energy, and mdi-icpton and £'di-icpton are the 
mass and the energy of the leptons originating from the Z decay. 
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Previous studies using simplified parametric detector simulations have shown the poten- 
tial of this technique [3]. 

Here, we study the prospects to measure the Higgs boson mass and cross section assuming 
its mass to be 120 GeV. At a ccntre-of-mass energy of 250 GeV, a full detector simulation of 
the process e"'"e~ ZH is performed using the Mokka software package, which simulates 
events in the LDC detector. The response from the sub-detectors is digitised as in a real 
experiment and is processed using the MarlinReco [4] reconstruction software. In addition, 
SM background processes are treated in the same way. 

The Z is reconstructed from its decays into electrons and muons. Algorithms are de- 
veloped to identify electrons and muons using the tracker and calorimeter information. A 
likelihood technique is used to separate signal events from the SM background processes 
with high efhciency. 



2 Experimental Conditions 

The study assumes a linear e~^e~ collider operating at a centre-of-mass energy of 250 GeV. 
This energy is chosen because the Higgs-strahlung cross section for the SM Higgs boson with 
a mass of 120 GeV reaches its maximal value. 

The statistics for signal and background events corresponds to a luminosity of 50fb~^. 
Signal events were generated using Pyhtia 6.4.11 [5]. Initial and final state bremsstrahlung 
and beamstrahlung are taken into account. To simulate bcamstrahlung the GuineaPig 1.4.1 
[5] program is used assuming ILC nominal beam parameters. Background events are pro- 
duced using in addition the event generators BHwide 1.04 [7j and Sherpa 1.0.10 [8] as 
listed in Table □ 

Signal and background events arc passed through a full detector simulation (Mokka) 
and are processed in the full reconstruction scheme of Marlin. A sketch of the simulation 
stages is shown in Figure [T] The LDC detector [9] is used for the simulation and reconstruc- 
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Figure 1: Scheme of simulation and analysis stages. Events of the different processes are generated 
using Pythia 6.4.11 , BHwide 1.04 and Sherpa 1.0.10. Beamstrahlung is treated by GuineaPig. 
Mokka and Marlin simulate and reconstruct events in the LDC detector. The analysis software 
is based on ROOT. 



tion. The vertex detector (VTX) consists of five layers of silicon pixel detectors. The main 
tracker is a TPC of about 3 m diameter and 4 m length supplemented by cylindrical silicon 
strip detectors (SIT) and forward strip and pixel detectors (FTD). The TPC is surrounded 
by the electromagnetic (ECAL) and hadronic (HCAL) calorimeter, which in turn are en- 
closed by the 4T magnet and the iron yoke. ECAL is a finely segmented silicon-tungsten 
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sandwich calorimeter. HCAL consists of a steel absorber structure with small scintillator 
tiles read out with silicon photomultipliers. When this analysis was performed no muon 
chamber was implemented in Mokka. Thus the separation of muons and pions in the par- 
ticle identification is done using the trackers and calorimeters only. The precision of the 
momentum measurement of the tracker system (TPC+VTX-I-SIT-I-FTD) is obtained to be 
Gpjpt = 7 ■ 10~^ • Pt [GeV] using the FuULDCTracking processor in MarlinReco. 



Process 


a[fb] 


iV(50 fb"') 


Generator 


1. 




HZ Xt 




15.0 


751 


Pythia 


2. 


e+e^ 






4144.5 


207223 


BHWIDE 


3. 


e+e^ 






4281.0 


214050 


Pythia 


4. 


e+e^ 






4182.0 


209100 


Pythia 


5. 


e+e^ — > 


w+w- 


Xe, Xfj., XepL 


5650.0 


282277 


Pythia 


6. 


e+e^ — > 


e+e-ff 




475.7 


23784 


Sherpa 


7. 


e'^e~ — ^■ 






359.4 


17970 


Sherpa 


8. 


e~^e~ —* 






24.6 


1231 


Sherpa 


9. 


e~^e^ 






7.2 


360 


Sherpa 


10. 








177.0 


8850 


Sherpa 



Table 1: The processes simulated for this study, their cross sections, the expected statistics for an 
integrated luminosity of 50 fb~^, and the generators used, i represents e, /i and / stands for r, v, q. 



3 Simulation Details and Analysis 

Signal and Background Processes 

The signal and background processes considered in the analysis, their cross sections and 
expected statistics at 50 fb~^, and the programs used to generate events are listed in Table[TJ 
For the generation of e+e~ events using BHwide, the following cuts are applied: the polar 
angle range is restricted to |cos?9iop| < 0.985, the electron energy is required to be Ee > 
10 GeV, the difference of the di-lepton mass and the Z mass is \mee — mzl < 40 GeV, and 
the recoil mass against the Z is in the range 90 GeV< mrccoii < 190 GeV. For event samples 
generated with Sherpa, cuts on the Icpton polar angle, | cosi9iop| < 0.985, the di-lcpton and 
di-parton mass, mee,qq > 10 GeV, and the energy of the final state fermions, -Efcrmion > 5 GeV, 
are applied. These cuts avoid divergences in the cross sections, reduce computing time, and 
would have no influence on the results. 

Analysis Strategy 

The signature of events from c+e^ ZH are two leptons of the same kind and opposite 
charge. The invariant mass of the two leptons must be in the vicinity of the Z mass. The 
dominant background is expected from the process e+e^ ZZ ^ £^£^X, which is simulated 
within the four-fermion processes 6 - 10 in Table [TJ Discriminating power is expected from 
the polar angle distribution of the Z. Processes 2 and 3 with two electrons or muons in the 
final state may be selected since initial state radiation leads to a radiative return and thus 
to an invariant mass near the Z. To distinguish them from the signal, the acoplanarity angle 
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can be used. The process 5 is a possible background due to its high cross section. The polar 
angle of the leptons can be used to distinguish the signal from this background. 

In the fist step of the analysis, we look for isolated electrons and muons in each event 
using a likelihood method. An electron is identified as an electromagnetic shower in the 
ECAL that matches the position predicted from a track as impact point into the ECAL. A 
muon is identified as a track matching to deposits in the ECAL and HCAL compatible with 
the expectations from a minimum ionising particle. 

For the second step, only events with a lepton pair of the same kind and with opposite 
charge are accepted. If there are several pairings possible, the one with an invariant mass 
closest to TOz is chosen. 

For further reduction of the background the following cuts are applied to: 

- the polar angles of the leptons: | cosz?iep| <0.95, 

- the difference between the invariant di-lcpton mass and mz: |TOdi-icpton ~ < 
30 GeV, 

- the lepton energy: i?icp>15 GeV, 

- the recoil mass: 90 GeV< mrccoii < 190 GeV, 

- the polar angle of the di-electron system: | cos-i^di-cicctronl < 0.90. 

The remaining events are analysed using likelihood density functions for the signal, the 
e"'"e~ — > 4f , and the e+c^ 2f background channels in the variables: acoplanarity angle of 
the two leptons, acolincarity angle of the two leptons, the di-lepton mass, the polar angle of 
the di-lepton system, the polar angles of the two leptons and the transverse momentum of 
the Z. 

For each event a likelihood is calculated characterising its compatibility with a signal 
event. A cut on this likelihood is set such that the quantity ^^g^ is minimised in the mass 
range from 119 GeV to 125 GeV. Here, S and B are the numbers of signal and background 
events, respectively, in the final sample. 

4 Results 

In the final sample, the signal selection efficiency is obtained to be 43.1% for the e+e" 
Zh e"'"e~X and 57.2% for the e+e" — > Zh — > ^"'"/j,~X channel, respectively. The recoil 
mass spectra are shown in Figure[2]for these processes. In both spectra a signal peak is seen 
on top of a moderate background. The signal has a smaller width in the di-muon channel, 
presumably because the muon track measurement is more precise due to less brcmsstrahlung 
in the material of the detectors. This is confirmed by the resolution function for the trans- 
verse momentum, which has a pronounced tail to lower reconstructed momenta for electrons. 



Cross Section Measurement 

The recoil mass spectra in Figure [2] are used to determine the cross sections for the processes 
e"'"e^ Zh e"*"e~X and e"'"e^ Zh —t /i+/i^X. The background originating from known 
SM processes is parametrised by a polynomial and kept constant in the fit to determine the 
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Figure 2: The recoil mass distributions of the selected e^e~ ZH events, left for Z — > e^e~ and 
right for Z — > fJ,'^ IJ.~ final states, respectively. The dark red distribution originates from the signal 
process and the light red distribution from the remaining background. 



amount of the signal. The signal is described using the following parametrisation, 

where mo is the central value of the peak. A a constant to describe the tail to larger values 
in the signal mass distribution and 1 — /? is the fraction of the tail. The tail to larger mass 
values in the signal is caused by bremsstrahlung and beamstrahlung. The former is well 
predicted from QED and the latter depends on the machine parameters. For known and 
reasonable stable machine parameters the shape and fraction of the tail can be determined, 
and we keep it constant in the fit, varying only toq and the number of events in the signal, 
-^signal- The cross section of the signal process is obtained from 

(T(process) = N^ignai / (Ce) , (3) 

where C is the integrated luminosity, and e is the signal selection efficiency. The results 
obtained are fT(e+e~ — > ZH) = 216.0 fb with an uncertainty of 20% using the di-electron 
final state and a(e"'"e~ ZH) = 219.7 fb with an uncertainty of 10% using the di-muon final 
state. Both results agree with the value of the cross section of 226.8 fb obtained from Pythia. 



An alternative method is to count all events, N, in the signal mass range from 119 GeV 
to 125 GeV and to subtract the background. The latter is obtained from a high statistics 
Monte Carlo simulation or from the integral over a parametrised background distribution in 
the same mass range, (B). The cross section is then given by 

(T(process) ^ {N - {B))/{Ce), (4) 

with an uncertainty of (± \/N / {N — (B) ) [%] . Using this method no assumption on the 
signal peak parametrisation is needed. The results obtained are compatible with the fit 
results given above. 
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Higgs boson mass 



To determine the Higgs boson mass from the spectra shown in Figure [5] a Ukehhood method 
is used. Several signal samples of high statistics with Higgs boson masses between 119 and 
121 GeV are generated and processed through the full simulation, reconstruction and analysis 
chain. The obtained spectra are parametrised using Formula ([2]). These parametrisations 
are then used in an unbinned likelihood fit to the simulated event samples shown in Figure 
[2] to determine the Higgs boson mass nih- The results are nih = 119.78 ± 0.42 GeV and 
ruh = 120.09 ± 0.12 GeV for the Z e+c^ and Z A^^M^ decays, respectively. 




Figure 3: Uncertainties on the Higgs boson recoil mass measurement (left) and on the Higgs boson 
production cross section (right) as a function of the centre-of-mass energy. A Monte Carlo toy 
model with parametrised momentum resolution is used. 



Estimate of the Optimal Centre-of-Mass Energy for the Measurements 

A Monte Carlo toy model is used to estimate the optimal centre-of-mass energy for the mea- 
surement of the Higgs boson mass and cross section. The resolution of the track momentum 
measurement is parametrised as (Jpjpt = 10^^ • pt [GeV]. A lepton identification and pair 
matching is performed as described above. For a Higgs boson mass of 120 GeV the recoil 
mass spectra are used to determine the Higgs boson mass and cross section for centre-of- 
mass energies between 210 GeV and 250 GeV. The estimated accuracies for the mass and 
cross section measurements arc shown in Figure [3] as a function of the centre-of-mass en- 
ergy assuming an integrated luminosity of 500 fb~^. The minimal uncertainty on the recoil 
mass occurs for i?cms = 220 GeV, in agreement with the results obtained in Ref . [TU] . The 
cross section uncertainty is minimal for E'cms = 240 GeV, becoming worse by about 20% at 
i;cms = 220GeV. 



5 Conclusions 



The recoil mass technique is a unique tool to determine the mass and cross section of the 
Higgs boson at the ILC. For the first time the prospects obtained from a full detector 
simulation and reconstruction are presented. Choosing the center-of-mass energy a few ten 
GeV above the kinematic threshold, here 250 GeV for ruh = 120 GeV, the Higgs boson mass 
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and cross section can be determined with an accuracy of 120 MeV and 9%, respectively, 
using only 50fb~^. To reach a similar accuracy at a center-of-mass energy of 350 GeV, an 
integrated himinosity 500 fb~^ is needed. 

The talk held at the LOWS is available under Ref. [TT] . 

References 

[1] S.L. Glashow, Nucl. Phys. 22 (1961) 579; S. Weinberg, Phys. Rev. Lett. 19 (1967) 1264; A. Salam, 
Elementary Particle Theory, edited by N. Svartholm (Almqvist and Wiksell, Stockholm, 1968), p. 367. 

[2] P.W. Higgs, Phys. Lett. 12 (1964) 132, Phys. Rev. Lett. 13 (1964) 508 and Phys. Rev. 145 (1966) 
1156 ; F. Englert and R. Brout, Phys. Rev. Lett. 13 (1964) 321; G.S. Guralnik, C.R. Hagen and 
T.W.B. Kibble, Phys. Rev. Lett. 13 (1964) 585. 

[3] P. Garcia- Abia and W. Lohmann, Eur. Phys. J. direct C 2 (2000) 2 [arXiv:hep-ex/9908065] . 

[4] ILC Software Portal: http://www-flc.desy.de/ilcsoft 

[5] T. Sjostrand, S. Mrenna and P. Skands, JHEP 0605 (2006) 026 [arXiv:hep-ph/0603175| . 

[6] D. Schulte, Universitat Hamburg (1996) 

[7] S. Jadach, W. Placzek and B. F. L. Ward, Phys. Lett. B 390 (1997) 298 [arXiv:hep-ph/9608412| . 

[8] T. G leisberg, S. Hoche, F . Krauss, A. Schalicke, S. Schumann and J. C. Winter, JHEP 0402 (2004) 
056 arXiv:hep-ph/0311263 . 

[9] DOD, LDC Working Group: http://www.ilcldc.org/documents/dod 

[10] F. Richard and P. Bambade, LAL-Orsay 07/03. 

[11] Slides: http : //ilcagenda. linearcollider . org/materialDisplay .py?contribId=159&materialId=sli 
desfeconf Id=1296 



LCWS/ILC2007 



